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This article briefly describes the nuance neutrino simulation software and outlines the program of the working 
group on neutrino event generators which met for the first time at the NUINT'01 meeting. 



1. HISTORY AND OVERVIEW 

The unrequited search for proton decay with 
large underground detectors, begun almost 20 
years ago, focused attention on describing at- 
mospheric neutrino reactions in detail, [jl] || Un- 
like typical accelerator beams of neutrinos, the 
spectrum of atmospheric neutrinos covers many 
decades of baseline and energy. From the lowest- 
energy contained interactions to the highest- 
energy upward-going muons, neutrino energies in 
a broad range (~100 MeV to >1 TeV) are sam- 
pled with roughly comparable rates. These fea- 
tures, which would eventually help reveal neu- 
trino mass and oscillation, demand a model of 
exclusive reactions valid for essentially any inter- 
acting neutrino energy. 

With the advent of Super-Kamiokande, the sta- 
tistical precision of atmospheric neutrino mea- 
surements increased dramatically and system- 
atic uncertainties (fluxes, cross-sections and cal- 
ibration) became increasingly important. The 
same situation is likely to arise in future long- 
baseline experiments with high-luminosity neu- 
trino beams. 

The nuance software package described here is 
the author's work-in-progress on the problem of 
modeling neutrino interactions]^] Although origi- 
nally used for atmospheric neutrino interactions, 
the program was designed to be useful in as many 
other applications as possible. 



1 This write-up reflects version 2.000 of the code. 



2. NEUTRINO PHYSICS MODEL 

Space does not permit a full elaboration of the 
cross-section calculations for each reaction chan- 
nel or graphical comparisons with experimental 
data, so a schematic outline with references must 
suffice. 

2.1. Generalities 

Nuance adopts a "divide and conquer" strat- 
egy. Models (with varying degrees of sophisti- 
cation) can be found in the literature for each 
general class of reaction. By summing the cross- 
sections and rates of all exclusive channels, and 
then adding (inclusive) deep-inelastic scattering 
inside appropriate kinematic limits, the total 
cross-sections and event rates are obtained. 

The units used by the program are MeV, grams 
(for densities) , centimeters and seconds. The user 
does not directly interact with cross-sections, but 
these are stored internally in units of 10 -36 cm 2 , 
or picobarns. 

Simulated elementary particles are indexed 
using PYTHIA[§ /Particle Data Group@] codes. 
Application-defined codes are used for nuclei in 
the few cases they are needed. Particle proper- 
ties and physical constants are generally taken 
from PYTHIA internal stuctures by default, al- 
though any parameter can be overridden at run- 
time. Where more recent values are available, 
they are taken from the Review of Particle Prop- 
erties. The PYTHIA particle stack is used to build 
the event in memory and process some decays. 

The Fermi gas model is used to simulate the 
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effect of a bound nucleon target (Fermi motion 
and Pauli blocking) . The implementation of this 
model is slightly different from process to process, 
however in general bound nucleons are given a 
uniform initial momentum density up to a (user- 
specified) maximum value, and a negative bind- 
ing energy. While different shells (with different 
momenta and binding energy) are supported, in 
practice the values Pf < 225 MeV/c and Et = 
—27 MeV (interpolated from a fit to electron- 
scattering dataQ) are used for all 16 nucleons 
in the canonical case of 16 O. A final state nu- 
cleon must exceed a threshold momentum value 
to exit the nucleus and allow the reaction to oc- 
cur. In the absence of final-state interaction, the 
(observable) outgoing nucleon momentum must 
therefore be greater than 225 MeV/c. Discussions 
at the Workshop made it abundantly clear that 
this hard cut-off is unphysical, however a more 
satisfactory presciption remains elusive. In most 
cases, nucleons below 225 MeV/c will not be de- 
tectable, but prudence is advised in interpreting 
low-energy recoil nucleons. 

Lepton masses are never neglected, and the ap- 
propriate vector-boson mass is always included in 
the propagator. 

2.2. Electron Scattering 

The only exactly-calculable (at first-order) neu- 
trino cross-section is unfortunately also the least 
important for neutrino energies above a few tens 
of MeV. For completeness, all purely leptonic neu- 
trino and anti-neutrino reactions are treated, in- 
cluding elastic scattering and inverse muon-decay 
(V e e~ — > /j,~v fJ/ ). The tree-level cross-section for- 
mulae are given in many particle physics text- 
books. || 

2.3. Quasi-elastic Scattering 

In nuance quasi-elastic scattering comprises 
both charged- and neutral-current two-body neu- 
trino reactions with nucleons. The relativistic 
Fermi gas model of Smith and Moniz provides the 
general framework for all such processes. [0 To 
ensure consistency between bound and free tar- 
gets, free nucleon cross-sections are also calcu- 
lated using the Smith-Moniz formalism, setting 
the binding energy equal to zero and expressing 



the initial nucleon momentum distribution as a 
delta function at zero. 

Identical form factors are used for both free and 
bound nucleons. The usual dipole parameteri- 
zation of vector and axial-vector form factors is 
adopted, with default values my = 0.840 GeV/c 2 
and vtia = 1-00 GeV/c 2 , respectively. For the 
induced pseudoscalar form factor, a parameter- 
ization calculated from lattice QCD|jll[] (which 
agrees extremely well with low-q 2 pion electro- 
production data) is used by default, although the 
simpler pion-pole expression from JTo) ] can be se- 
lected instead. Second-class currents are assumed 
to vanish, as required by the V — A electroweak 
theory and fundamental symmetries. 

The total and differential charged-current 
cross-sections of the Smith-Moniz model agree 
well with more sophisticated theoretical calcu- 
lations for neutrino energies above 50-100 MeV 
where continuum excitation of 16 O is dominant. 

Neutral-current two-body reactions with nucle- 
ons are usually invisible, but have been incor- 
porated into the framework of the Smith-Moniz 
model and included for completeness. Neutral- 
current nucleon form factors are specified by the 
electroweak theory and can be related to those 
for charged-current reactions. [fl2[ 

Similarly, the Smith-Moniz model has been 
extended to include charged-current, Cabibbo- 
suppressed hyperon production, following the 
treatment of Pais(l3) to account for the inelas- 
ticity of such reactions and the |A/| = \ rule. 

2.4. Resonant Processes 

For neutrino energies around 1 GeV and above, 
baryon resonances may be excited and subse- 
quently decay into a nucleon and one or more 
mesons. The most prominent resonance is the 
A(1232), but AT(1440) and a number of others can 
also contribute. Rein and Sehgal|Q have shown 
that these resonance-mediated channels can be 
described using harmonic oscillator quark wave- 
functions, symmetrized under SU(6) (extending 
the approach of Feynman, Kislinger and Ravn- 
dal|ll| to account for higher-mass resonant states 
and the isospin structure of the weak interac- 
tion). For hadronic masses above the A(1232), 
the data appear to require inclusion of interfer- 
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ence between resonances with identical isospin. 
Nuance adopts the Rein-Sehgal model, modified 
to account for improved knowledge of the mass 
spectrum since publication of their original pa- 
per, and incorporates non-strange resonances up 
to 2 GeV. Nucleon form factors for resonance pro- 
duction are assumed identical to those for quasi- 
elastic scattering, however Rein and Seghal ne- 
glect the pseudoscalar form factor of the nucleon, 
which (in principle) could be important for v T 
charged-current reactions.^ 

Rein and Sehgal consider only resonance de- 
cays into Ntt final states, but additional channels 
are possible. For completeness, the program at- 
tempts to model these by simply adding the ad- 
ditional decay modes so that the total branching 
fraction for each resonance sums to 100%. The 
procedure assumes the decay matrix element for 
other modes is identical to the Ntt case, which 
is probably wrong, but the contribution of these 
more exotic reactions to the total cross-section is 
extremely small. Unfortunately, due to the com- 
plexity of the Rein-Sehgal calculation, the large 
number of independent integration variables, and 
the increased number of channels, the bulk of the 
program's computing time is spent on these rela- 
tively unimportant reactions. 

For reactions on bound nucleons, the initial 
state is given a uniform Fermi momentum den- 
sity and a negative binding energy. Reactions 
in which resonance decay would result in a nu- 
cleon below the Fermi sea are Pauli-blocked and 
do not occur. Kim et al.]l6|, Singh and collab- 
oratorsJlT]] and Marteau et al.p^] have pointed 
out that the nuclear medium can modify the 
widths of resonances and allow the final-state re- 
action N*N — > NN. In-medium effects on the 
widths of resonances are not considered by the 
program, but the "pion-less A decay" reaction 
(a misnomer, since it seems relevant to J = I 
channels as well) can reduce the number of pions 
produced by 10-50% (the numbers suggested in 
the literature vary considerably). It is difficult to 
disentangle N*N — > NN from similar processes 
in which TV* — > Ntt decay occurs but the pion 

2 In practice, deep-inelastic scattering is the dominant v T 
charged-current reaction channel for most kinematically 
interesting neutrino energies. 



is subsequently absorbed, and it is not clear to 
what extent the former reaction is double-counted 
by final-state interactions. Nevertheless, Super- 
Kamiokande and K2K data appear to favor some 
suppression of resonant pion production. Absent 
a theoretically well-motivated and consistent de- 
scription of N* N — > NN in the Fermi gas model, 
by default nuance adopts an ad hoc 20% suppres- 
sion of pion production for 1% = ±i reactions and 
a 10% suppression for I 3 = ±§.[] Hopefully inter- 
action between the nuclear- and neutrino-physics 
communities fostered by the NUINT workshop 
will help clarify this important effect. 

2.5. Coherent and Diffractive Reactions 

In coherent reactions, neutrinos scatter from 
an entire nucleus rather than its individual con- 
stituents, with negligible energy transfer to the 
target. Coherent reactions typically produce a 
forward-going lepton and meson; both charged- 
and neutral-current reactions are possible, lead- 
ing to charged or neutral meson production, re- 
spectively. As a purely axial-vector process, the 
cross-sections for neutrinos and anti-neutrinos are 
equal, and the neutral-current cross-section is half 
as large as for charged currents. 

Nuance implements Rein and Sehgal's calcu- 
lation |nj of the cross-section for coherent pion- 
production in terms of the ttA forward scatter- 
ing cross-section. Rein and Sehgal's cross-section 
agrees with the handful of measurements avail- 
able, however the lowest-energy data available is 
at E v — 2 GeV, with an Al target. Using a 
far more detailed model of the nuclear physics, 
Kelkar et al.j2(| have stressed the importance of 
nuclear medium effects on the A which mediates 
these reactions, and predict a dramatic suppres- 
sion of the coherent cross-section for neutrino en- 
ergies around 1 GeV. Even with the larger cross- 
sections predicted by Rein and Sehgal, coherent 
pion production represents a small fraction of the 
incoherent, resonance- mediated pion-production 
cross-section for energies around 1 GeV, so the 
impact of this effect on the total pion production 
rate should be 10% or less. Data from the near 

3 Na'ively, a A ++ or A - can only participate in the 
N*N — > NN reaction with half as many nucleons as A+ 
and A . 



4 



detectors of K2K, and eventually MiniBooNE, 
should allow the coherent contribution to be sep- 
arated from the resonant channels based on the 
observed single-pion angular distribution. 

Diffractive reactions in nuance are analogous 
to coherent, except free protons rather than 16 O 
nuclei are the targets; the dynamics of the two 
reactions are identical. Rein's calculation pl| of 
the diffractive pion-production cross-section, us- 
ing a model similar to the coherent case, is im- 
plemented in nuance by the same routines which 
handle coherent reactions. The smaller size of 
the target reduces the cross-section correspond- 
ingly, by approximately a factor 16 a. Although 
the nuclear-medium effects cited in ^| do not 
apply to free nucleon targets, the cross-section 
for diffractive single-pion production is dwarfed 
by incoherent resonant channels and coherent re- 
actions with 16 O, hence it is included mainly for 
completeness. 

Coherent and diffractive production of vector 
mesons has been measured by several high-energy 
experiments, at the few per-mille level compared 
to the total charged-current rate. Cross-sections 
for these reactions, which are presently ignored in 
nuance, have been estimated based on the CVC 
hypothesis and vector meson dominance and ap- 
pear to describe the data reasonably well. For 
completeness, these channels will be added in a 
future version of the program. 

2.6. Deep-inelastic Scattering 

Deep-inelastic scattering is unique in that it is 
modeled as an inclusive, rather than exclusive, re- 
action between a neutrino and the parton consi- 
tituents of the nucleon. Elementary formula? are 
derived in most particle physics texts, however 
the usual expressions involve a number of unde- 
sirable kinematic approximations, including ne- 
glect of lepton (and target, in some cases) masses 
and the induced pseudoscalar form factor. The 
unpublished calculation of Roe[^3| avoids these 
limitations, and forms the kinematic basis for 
the treatment of this reaction in nuance. Addi- 
tional complications arise at low-energies, where 
the threshold for production of real heavy quarks 
(c,b) is important. The usual approach is to 
impose "slow rescaling"; here the calculation in 



nuance follows Leader and Predazzi.|24j 

Nuance interfaces to the PDFLIB package, |2^| 
allowing one of dozens of nucleon structure func- 
tion parameterizations to be selected at run-time. 
By default, the program uses the BEBCj26) op- 
tion, chosen because it was measured with neu- 
trinos and has a relatively low q 2 cut-off. Care is 
required in dealing with structure functions pro- 
vided by PDFLIB since requests for structure func- 
tions with \q 2 \ < \q 2 nin \ result in the values for 
q 2 = q^in being returned. In this case, the be- 
havior of the structure functions as q 2 — » is 
extrapolated according to the vector meson dom- 
inance model of PYTHIA.Q 

Shadowing, anti-shadowing, Fermi and "EMC" 
effects on the structure functions are neglected. 
Form factors are built from the structure func- 
tions assuming the Callen-Gross relation (F2 = 
2xFi) is exact and thereby also neglecting the 
longitudinal cross-section. 

Since the deep-inelastic cross-section includes 
resonant (and perhaps quasi-elastic) processes 
with low hadronic mass (which are already in- 
cluded via exclusive channels) the kinematic 
limits of integration must be chosen to avoid 
double-counting. As extensively discussed at 
NUINT'01, in real-life there is no distinct cut- 
off between "resonant" and deep-inelastic scat- 
tering, but rather a smooth transition from scat- 
tering off hadrons to scattering off quarks, as the 
vector boson probe begins to resolve the internal 
structure of the target. The present version of 
the nuance attempts to square this circle by inte- 
grating the deep-inelastic cross-section over the 
limits \q 2 \ > 1 (GeV/c) 2 .OR. W > 2 GeVfl 
These values are chosen to make the differential 
cross-section roughly continuous across the artifi- 
cial boundary at W = 2 GeV, the upper limit for 
resonant reactions in the program, although it is 
only partially successful. 

Because the deep-inelastic scattering cross- 
section is inclusive, these channels are different 
from others in the program, where the final-state 
particles are uniquely determined by the reaction 
itself. For deep-inelastic scattering, the interact- 



4 The additional constraint W > rrajv+m-ir is also imposed, 
regardless of q 2 . 
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ing quark is determined by a weighting scheme 
from the LEPTO program^] and transformed into 
a different flavor (for charged-current reactions). 
Parton showers and fragmentation of the outgo- 
ing quark are treated by a version of LEPTO which 
has been extensively modified to ensure conser- 
vation of energy and momentum and has kine- 
matic cuts adjusted so |<?^ in = 1CP 3 (GeV/c) 2 
and W^ in = 2 GeV 2 . LEPTO itself uses PYTHIA 
for some of its work. 

Because low-mass and small- q 2 reactions are 
sometimes generated, the LUND string fragmen- 
tation model occasionally fails to produce a final 
state. In this case KNO scaling Q is used to de- 
termine a final-state using measured multiplicity 
distributions and fragmentation functions. p9| 

Deep-inelastic scattering is undoubtedly the 
regime where nuance currently lags behind the 
state-of-the-art in other generators dedicated to 
high-energy neutrino physics. The NUINT meet- 
ing has already borne fruit by bringing the issue of 
quark/hadron duality to the forefront and elicit- 
ing at least one clear prescription for how to rem- 
edy these shortcomings. This approach will be 
implemented in a future version of the program. 

2.7. Nuclear Processes 

The program uses a model of final state inter- 
action in the nucleus originally developed for the 
1MB experiment. pH At present the model is spe- 
cific to 16 O, although it could be adapted to other 
nuclei. Primary interactions are assigned a start- 
ing position in the nucleus according to the mea- 
sured density distribution for 16 0.j3^] Hadrons 
are then tracked through the nucleus in 0.2 fm 
steps, treating the nucleus as an isoscalar sphere 
of nuclear matter with radially-dependent den- 
sity and Fermi momentum. Single-nucleon cross- 
sections and local density are used to calculate the 
interaction probability during each step. Interac- 
tions resulting in a nucleon with momentum be- 
low the Fermi sea are Pauli-blocked and ignored. 

Measured cross-sections and angular distribu- 
tions are used for ir — N and N — N reactions. |p3[ 
Angular distributions for elastic reactions are cal- 
culated from a global phase-shift analysis of world 
data.|34| Inelastic reactions involving up to five 
particles are possible. For three-body ir — N reac- 



tions, a AI = | dominance model is adopted: |35f 
in other simple resonance model is as- 

sumed. For reactions without available data, 
cross-sections are inferred from isospin symmetry. 

Interactions of kaons with kinetic energy up to 
1 GeV are also simulated. Cross-sections and an- 
gular distributions for all two-body K — N re- 
actions are calculated using a partial-wave anal- 
ysis. |Q Neutral kaons are treated as 50% Ks 
(which decay immediately) and 50% Kl- Kl 
cross-sections are calculated as an equal mixture 
of if and K , and any interaction has a 50% 
chance of regenerating a Ks which (again) im- 
mediately decays. 

Hyperon interactions inside the nucleus are not 
simulated and p mesons decay rapidly before any 
interaction is possible. Elastic scattering and 
pion-production cross-sections J|lj compete with 
decay for u> and rj mesons inside the nucleus. 

The 1MB cascade model has been tested by 
running the algorithm in a mode which simulates 
7T, p, K — 16 scattering by starting a hadron out- 
side the nucleus and stepping it through. The 
pion absorption cross-section is tuned to repro- 
duce measurements; for other channels, the re- 
sults of the simulation reproduce scattering data 
for hadrons on 16 and 12 C extremely well for ki- 
netic energies up to 2 GeV. For pions, the model 
also agrees with the less detailed (but more ele- 
gant) analytical approach of Adler, Nussinov and 
Paschos.|| 

One flaw in such superposition models was 
pointed out at NUINT '01, namely the failure to 
account for "formation zones" , which suppress 
the final-state interactions of hadrons with ener- 
gies greater than a few GeV. (3^] In the present 
version of the program, cross-sections are as- 
sumed to be constant for energies above 2 GeV. 

After ejection of a nucleon recoiling from a neu- 
trino interaction, de-excitation and/or break-up 
of the residual nucleus can result in emission of 
one or two few-MeV 7 along with possible evap- 
oration of low-energy nucleons. |4(| While the de- 
excitation products are usually quite low in en- 
ergy, their signature in a water detector is suffi- 
cient to shift the average reconstructed mass of tt° 
produced in neutral current reactions by as much 
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as 5 MeV. The model of de-excitation currently 
used is specific to 16 O, and should be generalized. 

Once all neutrino interaction products have es- 
caped the nucleus or been absorbed, and the nu- 
cleus has de-excited, nuance completes its simu- 
lation of the event by calling PYTHIA to process 
decays of short-lived particles outside the nucleus. 

2.8. Utilities 

Nuance provides additional services for spe- 
cialized applications, and uses others which are 
publicly available. The program includes a 
package for modulating input fluxes under 3- 
component neutrino oscillations in matter with 
CP-violation |4l| (including a realistic "onion- 
skin" model of the Earth's density profile). In 
practice, one often prefers to generate unoscil- 
lated (charged- plus neutral-current) and "fully- 
oscillated" (charged-current only; P osc = 1, in- 
dependent of energy) samples, and then reweight 
them according to different hypotheses after the 
fact. The program will also generate nucleon de- 
cays in a user-specified channel, using the same 
bound nucleon and final-state interaction model 
applied to neutrino interactions. 

Decays of r leptons are handled by the TAUDLA 
package. p2[ The polarization of primary r lep- 
tons produced by neutrinos is calculated using the 
form factors for quasi-elastic and deep-inelastic 
reactions. p3[ In resonant and diffractive pro- 
cesses, the appropriate form factors are unclear 
and/or difficult to calculate, so the tau is assumed 
to be completely polarized. 

The program can also simulate entering 
neutrino-induced muons produced in the mate- 
rial outside a detector. For transport of energetic 
muons to the detector surface, the PRDPMU pack- 
age[|4| is used. Note that small-angle Coulomb 
scattering of muons in flight is neglected for com- 
putational efficiency. 

A stand-alone utility program (nuplot) is also 
provided to translate the calculated cross-sections 
and rates from the internal format of a sequen- 
tial access ZEBRA/FZ file created by the program 
into an Ntuple for inspection. Using this util- 
ity, energy-dependent cross-sections for any of the 
many exclusive reactions can be individually plot- 
ted. 



2.9. Validity and Limitations 

Unphysical final-state nucleon energies in the 
Fermi gas model, the 16 0-specific description of 
final-state interactions and a number of other 
problems and simplifications have been discussed 
above. 

Because the program is written in F0RTRAN77 
and uses ZEBRA |45|] for dynamical memory man- 
agement, single precision calculations are per- 
formed in most cases (PYTHIA uses double pre- 
cision calculations internally). This is an unsat- 
isfactory solution, since it requires considerable 
care in coding to preserve numerical precision 
(such problems typically appear for neutrino en- 
ergies approaching 1 TeV) . A future version of the 
program, using a language with native dynamic 
allocation (such as FDRTRAN90 or C++), will even- 
tually remove this handicap. 

3. INTERFACE 

The program was intended to be flexible and 
easy to use, although these requirements are fre- 
quently at odds, if not mutually exclusive. To 
allow wider access to the code, it has been posted 
online at http://nuint.ps.uci.edu/nuance. 
The code is developed and maintained under the 
Windows operating system, but is fully compat- 
ible with others. The CMT configuration man- 
agement tool|46) provides a platform-independent 
mechanism to build the program and its con- 
stituent libraries automatically. The main pro- 
gram itself consists of only eight lines of exe- 
cutable code, hence the libraries that do the ac- 
tual work can be easily slaved to an alternative 
user-written program. 

The program is steered using a ZEBRA/TzQ 
(text) cards file, although frequently-changed op- 
tions can also be specified on the command line. 
Default physical parameters are stored in a sec- 
ond cards file named nuance_def aults . cards. 
Values specified in the job-specific file provided by 
the user override the defaults. The user can also 
modify any PYTHIA parameter via a data card. 

3.1. Setup 

The user must describe the geometry and com- 
position of the target before running the program. 
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The description syntax is reminiscent of GEANT, 
but somewhat simpler. A target consists of a one 
or more volumes (sphere, cylinder or box) which 
can be nested, but presently must all be centered 
at the origin. Each volume contains a mixture of 
one or more materials, with user-specified density. 
Materials are in turn composed of atoms, which 
are built from neutrons and protons in shells with 
specific Fermi momenta and binding energy. Each 
atom is assumed to include a number of electrons 
equal to the number of protons. 

3.2. Fluxes 

The program can generate events for a mono- 
energetic, single-flavor beam of neutrinos, how- 
ever for most applications the user must also pro- 
vide a flux of neutrinos. Drivers compatible with 
several atmospheric flux calculation formats are 
included. 

An accelerator beam can be described via an 
HBOOK file containing one or more histograms 
(each associated with a particular neutrino fla- 
vor) in units of neutrinos per bin per cm 2 per unit 
luminosity. These histograms must have iden- 
tical upper and lower limits, but the user may 
specify whether the limits are interpreted as MeV 
or GeV. Two scaling factors, one linear and one 
quadratic, allow the luminosity units and base- 
line of the beam (respectively) to be changed at 
run-time. The beam direction in local coordi- 
nates is also adjustable. The user may apply 
multi-quadric smoothing on the input histogram 
to eliminate binning effects, or directly generate 
neutrinos from the raw values. 

Finally, a simple model of thermal neutrino 
emission (for supernova?) [^7j is built into pro- 
gram, allowing the user to specify a total lumi- 
nosity and cooling time for one or more neutrino 
flavors. In this special case, the fluxes are time- 
dependent and interaction rates are recalculated 
as the simulation progresses and the source cools. 

3.3. Cross-sections and Rates 

Having described what is to simulated, the user 
must calculate neutrino cross-sections and inter- 
action rates. The caculation is time-consuming (a 
significant fraction of a day on a 2 GHz Pentium 
machine) due to the multi-dimensional numerical 



integrations involved, but need only be carried 
out once for a particular energy range. The cross- 
sections and rates are stored in a ZEBRA/FZ file for 
use during event generation; the interaction rates 
can be recalculated rapidly for a different target 
size or different flux with identical energy lim- 
its by reusing the stored cross-sections. Various 
classes of reactions, or specific neutrino flavors, 
may be selectively processed or ignored via flags 
given on the command line. 

Multi-dimensional integration of the rate and 
cross-section for each channel is performed by 
the Monte Carlo routine MISER or by cas- 
caded Romberg integration. [Q In both cases, a 
user-selectable tolerance on the estimated frac- 
tional error of the integral (0.3% by default) is 
checked to determine whether convergence has 
been achieved. Some kinematic integration vari- 
ables are transformed internally to accelerate con- 
vergence. The integral over neutrino energy nec- 
essary to compute total interaction rates can be 
done linearly or logarithmically. To ensure con- 
sistency between calculated rates and simulated 
data, the same code is used for both cross-section 
calculations and event generation. 

3.4. Event Generation 

Once cross-sections and rates are calculated 
and stored, event generation proceeds very 
rapidly; when generated event vectors are writ- 
ten to disk, the program is I/O-limited. The user 
indicates the exposure time in units of seconds 
or years (for atmospheric and supernova fluxes, 
respectively) or units of luminosity (for neutrino 
beams). A fixed number of events to generate 
may alternatively be specified. As in rate cal- 
culation, selected classes of reactions or neutrino 
flavors may be specified or suppressed. 

Each reaction channel is tracked separately to 
determine how much simulation time elapses be- 
tween events. After determining which channel 
produces the next event, a neutrino energy is gen- 
erated from the cumulative differential rate distri- 
bution for the reaction stored on the FZ file. A set 
of transformed kinematic variables is created ran- 
domly within the allowed limits and the differen- 
tial cross-section is calculated. The trial configu- 
ration is accepted or rejected by comparing a ran- 



8 



dom number to the ratio of the differential cross- 
section to the maximum differential cross-section 
for the same neutrino energy encountered during 
the rate integration phase (a table of a max vs. E v 
for each channel is also stored in the FZ file). If 
a particular set of kinematic variables is rejected, 
a new set is generated for the same reaction and 
neutrino energy until an accepted configuration 
is obtained. Using an accepted set of kinematic 
variables, four-vectors for the reaction's primary 
outgoing particles are generated. These particles 
are stepped through the target nucleus and/or de- 
cayed, if appropriate, and the final outgoing event 
vectors are written to disk. 

For obvious reasons, the program is a voracious 
consumer of random numbers. The RANLUX pack- 
age [[l9] is used everywhere (including external li- 
braries, which are modified to call it instead of 
their native random number generators). The 
user should provide a 32-bit integer seed in the 
steering cards or on the command line to ensure 
that different jobs produce unique results. 

The program writes its results in either (or 
both) of two formats as requested in the steering 
cards or on the command line: a human-readable 
text file or an HBOOK N-tuple. The text file lists 
not only the initial and final-state tracks, but also 
the primary outgoing particles prior to nuclear 
interactions and decays. The reaction channel, 
elapsed simulation time, neutrino flux at the gen- 
erated angle/energy and various internal or infor- 
mational kinematic variables (N-tuple only) are 
also recorded. Total charged- and neutral-current 
cross-sections for each neutrino species and tar- 
get combination are written to the HBOOK file in 
a separate N-tuple. 

4. TOWARD A UNIVERSAL MODEL 

4.1. Motivation 

As outlined in the preceding sections, accu- 
rately modeling neutrino interactions presents a 
unique challenge, combining the complexity of the 
nuclear many-body problem with the formidable 
uncertainties and sparse experimental guidance 
typical of neutrino physics. At the same time, 
discovery of neutrino oscillation has opened a new 
and inviting window on physics beyond the Stan- 



dard Model, which ambitious future long-baseline 
experiments now being planned aim to exploit. 

To date, each new experiment has had to solve 
the problem of simulating neutrino interactions 
independently. The literature contains a patch- 
work of disconnected and sometimes incompat- 
ible models. Piecing them together involves an 
equal mixture of intuition, folklore, and guess- 
work; presently, everyone working on the prob- 
lem must make their own guesses in a vacuum, 
largely by trial and error. Without a common 
point of reference, there is no easy way to test the 
validity of one's assumptions and no avenue for 
one group's progress or ideas to reach others. In 
short, there is virtually unlimited room for consol- 
idation and incremental improvement in the state 
of the art, but no mechanism for it to occur. One 
can only imagine the chaos which would prevail 
in pp or e + e~ physics if the LUND parton shower 
and string fragmentation code were not univer- 
sally available and each experiment were instead 
forced to recreate something similar on their own, 
yet this is precisely the present situation in neu- 
trino physics. 

Further, as interest in neutrino physics contin- 
ues to grow, more people enter the field and more 
future experiments are explored, the need for a 
general model only increases. On the other hand, 
the investment of time required to create a reli- 
able simulation is prohibitive for those interested 
only in exploring future possibilities rather than 
operating an approved experiment. Most "pro- 
prietary" code developed by a particular experi- 
mental collaboration tends to be specific to their 
own running conditions, and many groups are re- 
luctant to release internal software. 

4.2. The NUINT Working Group 

The first NUINT meeting was a watershed in 
the history of neutrino physics. In addition to 
bringing together a diverse group from across 
the spectrum of nuclear and particle physics to 
offer their specialized insights and expertise, it 
also united for the first time a large subset of 
the world's neutrino simulation experts, many of 
whom have labored in isolation for years. For me, 
this opportunity to compare notes and exchange 
ideas was the high-point of an already- fascinating 
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workshop. 

Agreement on the need for more comprehen- 
sive and widely available software tools was 
unanimously expressed, along with a willing- 
ness to work together toward them in the fu- 
ture. The first step will be to compare our re- 
sults in a number of agreed "benchmark" cases. 
This work will highlight the most important ar- 
eas of ambiguity and uncertainty, and hopefully 
motivate renewed theoretical attention on long- 
neglected but essential points. A central web site 
(http : //nuint . ps . uci . edu) linked to the home- 
pages of each participant will allow the working 
group and the community as a whole to track the 
progress of this effort and provide feedback. 

For the longer term, the efforts of many will be 
required to produce a carefully-tested and uni- 
versal model of neutrino interactions. In addi- 
tion to purely technical considerations, theoreti- 
cal guidance and new experimental data will be 
vital. Still, with the success of NUINT'01 and 
the promise of renewed and expanded collabora- 
tion punctuated and reinforced by future NUINT 
workshops, it is not too optimistic to hope that 
within a relatively few years, members of the neu- 
trino physics community will finally have at their 
disposal a software tool equal to the challenges 
and opportunities they face. 
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